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Recent Developments in Quantum Vacuum Energy for Confined Fields 
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Quantum vacuum energy entered hadronic physics through the zero-point energy pa- 
rameter introduced into the bag model. Estimates of this parameter led to apparent 
discordance with phenomenological fits. More serious were divergences which were omit- 
ted in an ad hoc manner. New developments in understanding Casimir self-stresses, and 
the nature of surface divergences, promise to render the situation clearcut. 
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1. Introduction 

In the bag modelPthe normal vacuum is a perfect color magnetic conductor, that 
is, the color magnetic permeability fj, is infinite, while the vacuum in the interior of 
the bag is characterized by ji = 1 . This implies that the color electric and magnetic 
fields are confined to the interior of the bag, and that they satisfy the following 
boundary conditions on its surface S: n • E|s = 0, n x B|s = 0, where n is a unit 
normal to S. Now, even in an "empty" bag (i.e., one containing no quarks) there 
will be nonzero fields present because of quantum fluctuations. This gives rise to a 
zero-point or Casimir energy. 
The bag-model Lagrangian is 

£ba g = (£ D + £ym - B)r](ipip), (1) 

where £d,ym are the Dirac and Yang-Mills Lagrangians, respectively, B is the bag 
constant, and rj is the unit step function. Variation of £bag leads to the following 
linear and quadratic boundary conditions (n M is the outwardly directed normal) 
—ijntp = tp and n a G a P = 0, and — jG 2 — ^(nd)ipip — B, in terms of its quark 
contribution. 

The MIT bag model consists in solving the Dirac equation subject to the bound- 
ary conditions, which in the simplest realization refer to a spherical cavity approx- 
imation. The interactions, through the gauge fields, are treated perturbatively in 
the approximation that the strong coupling is regarded as small, a s -C 1. The radial 
wavefunctions are then expressed in terms of spherical Bessel functions. The results 
of fits, shown in Table 1 are essentially identical with those found by DeGrand et 
al.pl Model e. The following comments can be offered concerning these fits. 
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Table 1. Fits to hadron masses for various bag models. 



Parameters 


Experimental 


Model 


Model 


Model 


and Masses 


Values 


A 


D 


E 


Z 




-1.0 


1.0 


1.8327 


B 1 / 4 (GeV) 
F (GeV 2 ) 




0.19526 


0.16515 


0.14535 




-0.21923 


-0.064729 





OLs 




2.16066 


2.17932 


2.19897 


m a (GeV) 




0.2807 


0.2802 


0.2797 


R N (GeV" 1 ) 




5.06 


5.03 


5.01 


8m (GeV) 




0.1325 


0.1343 


0.1408 


niK 


0.138 


0.3610 


0.3229 


0.2833 


rn n 


0.549 


0.645 


0.6105 


0.5583 


m v , 


0.958 


0.5031 


0.4702 


0.4209 


m K 


0.4957 


0.5827 


0.5489 


0.4982 


m N 


0.9389 


0.9389 


0.9389 


0.9389 


rn A 


1.1156 


1.1022 


1.1017 


1.1012 


m E 


1.1931 


1.1426 


1.1426 


1.1422 


mw 


1.13181 


1.2885 


1.2876 


1.1867 


m p 


0.769 


0.7826 


0.7826 


0.7826 


mu, 


0.7826 


0.7826 


0.7826 


0.7826 




1.0195 


1.0733 


1.0711 


1.0688 


m K * 


0.8921 


0.9251 


0.9246 


0.9229 


m A 


1.232 


1.232 


1.232 


1.232 


m s » 


1.3839 


1.3760 


1.3761 


1.3762 


m=t 


1.5334 


1.5229 


1.5229 


1.5230 




1.6725 


1.6725 


1.6725 


1.6725 



• The Z = — 1 (repulsive) estimate (my best guess in 1983*1), model A, agrees about 
as well with the data as does the optimal bag model E. 

• All the models agree more closely with each other than with the data. 

• Chiral bag models can be used to improve the fits with the pseudoscalars. 

• For models A and D a constant force term is included: (the bag radius is a) 

H = BV + crA + Fa+ - + H qg . (2) 
a 

• Some details of how Z is calculated will follow. 



1.1. Gluon and quark condensates 

Due to zero-point fluctuations, confinement will result in a nonzero expectation 
value for the squared field strength, (S — > is a time-splitting regulator) 

(G\r))=' ^±(21 + 1)2 r^ ^l^ a \ (3) 



Similarly, we obtain the following expression for the quark condensate 

-sf(xr/a) + sf +1 (xr/a) 

4{i -f l) j ax a 

1=0 



mr)) = -j^fxi + 1) r (4 , 

4?r ar J-°° sf(x)+sf +1 (x) 
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Fig. 1. Magnitude of the quark and gluon condensates as functions of the distance from the 
center of the bag. The values shown are for a single color, and, for the quarks, for a single flavor 
and helicity state. 

Here s; and e/ are modified spherical Bessel functions defined in (|1U|) below. Details 
of these calculations were given in Ref. 3 Graphs of these condensate densities 
are shown in Fig. ^ The central values seem to agree with the observed values of 
the quark and gluon condensates^ but they diverge as the surface of the bag is 
approached. 

2. Casimir Forces on Spheres 

The calculations presented here were carried out in response to the recent program 

of the MIT group. They rediscovered irremovable divergences in the Casimir energy 

|7T| 

for a circle first discovered by Sen in 1981," but then found divergences in the case 
of a spherical surface, thereby casting doubt on the validity of the Boyer calculation 
for the Casimir energy of a perfectly conducting spherical shell0 Some of their 
results, as we shall see, are spurious, and the rest had been earlier discovered by 
the Leipzig groupie However, the MIT group's work has been valuable in sparking 
new investigations of the problems of surface energies and divergences.^ 

2.1. S -function Potential 

We now carry out a calculation^^ of the zero-point energy in three spatial dimen- 
sions, with a radially symmetric singular background 

C int ^- l -- a 5{r~a)<t?(x) : (5) 
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which would correspond to a Dirichlet shell in the limit A — > oo. The time- Fourier 
transformed Green's function satisfies the equation (k 2 = — uj 2 ) 



A 



-V 2 + n 2 + -S(r - a) 
a 



G(r,r') =<5(r-r') 
We write G in terms of a reduced Green's function 



(G) 



(7) 



where gi satisfies 



1 d 2 d 
r 2 dr dr 



1(1 + 1) 2 A., 
V „ ; + k 2 + -S(r - a) 



gi(r,r') = -^S(r - r'). 



The solution inside the sphere (0 < r,r' < a) is 
1 



A 2 si{Kr)si(Kr) 

ez(Kr>)sz(/tr<) e, («o)— — — - — — — - 

Ka 1 + -^si(Ka)ei{Ka) 



Here we have introduced the modified Riccati-Bessel functions, 

F2x 



si(x) 



h+i/2 {%), ei(x) 



Ki +1/2 {x) 



(8) 



(9) 



(10) 



Note that this reduces to the expected result, vanishing as r — > a, in the limit of 
strong coupling. When both points are outside the sphere, r, r 1 > a, we obtain a 
similar result, with e; <-* s; in the second term. 



2.2. Pressure on sphere 

Now we want to get the radial-radial component of the stress tensor to extract the 
pressure on the sphere, which is obtained by applying the operator 



d r d r > - i(-<9°d /0 



d r d r i 



K 2 - 



1(1+1) 



(11) 



to the Green's function, where in the last term we have averaged over the surface 
of the sphere. In this way we find, from the discontinuity of (T rr ) across the r = a 
surface, the net stress, which can alternatively be obtained by differentiating with 
respect to a, with A/a fixed, the total Casimir energy 



1 °° f°° d 

E = -— J" (21 + 1) / dxx — \n[l + \I v (x)K v (x)] 
2va Jo dx 



(12) 



This expression, upon integration by parts, coincides with that given recently by 
Bartonpl and was first analyzed in detail by ScandurraU^ 
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2.3. Strong coupling 

For strong coupling, the above energy reduces to the well-known expression for the 
Casimir energy of a massless scalar field inside and outside a sphere upon which 
Dirichlet boundary conditions are imposed, that is, that the field must vanish at 
r = a: [y = I + 1/2) 



1 °° f 00 d 

lim E=-- — V(2/ + l) / dxx — \n\lJx)K v 

i=o JV 



(x)} , (13) 



because multiplying the argument of the logarithm by a power of x is with out ef- 
fect, corresponding to a contact term. This may be evaluated numerically!^ E TE = 
0.002817/a, which is much smaller than the Boyer result for electrodynamics 
e em _ o.04618/a, although both are repulsive. 

2.4. Weak coupling 

The opposite limit is of interest here. The expansion of the logarithm in (|12(l is 
immediate for small A. The first term, of order A, is evidently divergent, but ir- 
relevant, since that may be removed by renormalization of the tadpole graph. In 
contradistinction to the claim of of the MIT group^the order A 2 term is finite. That 
term is 

\2 °° roc j 

S(A2) = i^E^ + 1 ) / dxx-[I l+1/2 {x)K l+1/2 (x)f. (14) 
;=o 

The sum on I can be carried out using a trick due to KlichE] The sum rule 

oo 

V(2Z + l)ei{x)si{y)P l {co^6) = ^V", (15) 

P 

1=0 F 



where p — \/ x 2 + y 2 — 2xy cos 8, is squared, and then integrated over 9. In this way 
we learn that 

^2(2l + l)e"(x) S f(x) = ^- ^e-. (16) 

U 2 J ° w 

Although this integral is divergent, because we did not integrate (|12fl by parts, that 
divergence does not contribute: 

#)4r^4r-e-"=-, (17) 

47ra Jo 2 dx J w 32™ 

which is exactly the result I had found earlier,^! based on the following formula for 
a hypersphcrc in D dimensions: 

p(a 2 ) A^ r(^i)r(j-3/2)r(i-D/2) 

na 2 1 + 2D [T{D/2)] 2 ' 1 ' 

which exhibits poles when D is even, where the Casimir energy is known to 
diverge.E3 
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2.5. Divergences 



However, before we wax too euphoric, we recognize that the order A 3 term appears 
logarithmically divergent, just as the MIT group claimedPlTo study the behavior of 
the sum for large values of I, we can use the uniform asymptotic expansion (Debye 
expansion), 



I v {x)K v {x)~ — 



1 A(t) B(t) 



(19) 



Here x — vz, and t = 1/y/l + Z 2 . The functions A and B, etc., are polynomials in 
t. We now insert this into the energy expression and expand not in A but in v; the 
leading term is 

(A3) A 3 f,l dz A 3 

24™ j^vj (1 + ^ 2 ) 3 / 2 ~ 247ro U j ' ( ' 

Although the frequency integral is finite, the angular momentum sum is divergent. 
The appearance here of the divergent £(1) seems to signal an insuperable barrier 
to extraction of a finite Casimir energy for finite A. The situation is different in the 
limit A — > oo. 

This divergence has been known for many years, and was first calculated explic- 
itly in 1998 by Bordag et al.P where the second heat kernel coefficient gave 

E~-£-K s^O. (21) 
48™ s 

A possible way of dealing with this divergence was advocated in Scandurra.^l Very 
recently, Bordag and VassilevicfP^ have reanalyzed such problems from the heat 
kernel approach. They show that this 0(A 3 ) divergence corresponds to a surface 
tension counterterm, an idea proposed by m in 1980 in connection with the 
zero-point energy contribution to the bag model. Such a surface term corresponds 
to A/a fix ed, which then necessarily implies a divergence of order A 3 . Bordag and 
Vassilevicli^' argue that it is perfectly appropriate to render this divergence finite 
by renormalization. 



2.6. Boundary layer energy 

Here we 

showPl 

that the surface energy can be interpreted as the bulk energy of 
the boundary layer. We do this by considering a scalar field in 1 + 1 + d dimensions 
interacting with the background £i n t = — ■|0 2 cr, where 

h, 
0, 



a(x) 



2 < X < 2 , 



with the property that hS 
dx 2 



otherwise, 

1. The reduced Green's function satisfies (k 2 



\o{x) 



g(x, x') = S(x — x'). 



(22) 

kl-LU 2 ) 

(23) 
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This may be easily solved in the region of the slab, — | < x < |, («/ = \J n 2 + Xh) 
1 



g{x,x) = 



2k 1 



,-k'\x-x'\ , J_ 

A 



(n 12 - k 2 ) cosh /c'(x + x') 
+ (k' - n) 2 e~ K ' s cosh k' (a; - x') 
A = 2kk' cosh k' 5+ (k 2 + k' 2 ) sum k'<5. 



(24) 



(25) 

We generalize the considerations of Graham and OlumpH by considering the 
stress tensor with an arbitrary conformal term, 



= d^Wt - ^g^(d x <j)d x (j) + Xhcf 2 ) - a{d»d v - g^d 2 )<j) 2 . 
We get the following form for the energy density within the slab, 

2 -d-2 7r -(dfl)/2 f oo dKK d 



(26) 



T 



oo 



r((d + 3)/2) J K 'A 
- («' - K) 2 e 



{k' 2 - k 2 ) [(1 - 4a) (1 + d)n' 2 - k 2 ] cosh2 K 'a; 

(27) 



From this we can calculate the behavior of the energy density as the boundary is 
approached from the inside: 



T 



■oo 



T(d+l)Xh 



l-4a(d+l)/d 



\x\ -» 5/2. 



(28) 



2 d+4 7r (d+i)/2 r (( d + 3 )/2) (5-2\x\) d 

For d — 2 for example, this agrees with the result found by Graham and Olum for 
a = 0: 

Aft (l-6a) , 6 



00 



(29) 



96tt 2 (<5/2- |.T|) rf ' 1 

Note that, as we expect, this surface divergence vanishes for the conformal stress 
tensor, where a — d/4(d + 1). (There will be subleading divergences if d > 2.) 

We can also calculate the energy density on the other side of the boundary, from 
the Green's function for x,x' < —5/2, 



g(x,x') = — 

2k 



& -k\x-x'\ _ e K(i+i'+S)^/2 _ 



sinh k'5 



and the corresponding energy density is given by 



T 00 = 



d(l-4a(d + l)/d) 

" 2 d+2 7r (d+l)/2 r (( d+ 3)/2) J Q 



dn k 



d+l i K ' K ) „1k(x+&/2) 



A 



(30) 



sinh k 5, (31) 



which vanishes if the conformal value of a is used. The divergent term, as x — » —5/2, 
is just the negative of that found on the inside. This is why, when the total energy is 
computed by integrating the energy density, it is finite for d < 2, and independent 
of a. The divergence encountered for d — 2 may be handled by renormalization of 
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the interaction potential. In the limit as h — > oo, hS = 1, we recover the divergent 
expression for a single interface 

1 f°° A 

^oo Es = 2*+W+WT((d + 3)/2) Jo ^ (32) 
Therefore, surface divergences have an illusory character. 



2.7. Dielectric Spheres 

The Casimir self-stress on a uniform dielectric sphere was first worked out by me 
m 1979.H1 It was generalized to the case when both electric permittivity and mag- 
netic permeability are present in 1997.2^ The result for the pressure (x — ^/s]I\y\, 
i 



x 



/e'fi'\y\ where e',fj/ are the interior, and s,fi are the exterior, values of the 



permittivity and the permeability) is 

P =-^r d ^± 2 -i±i{x± M 

2a 4 J /_ 00 2 7 r ^ in \ dx 

+ ^ x '[s'i{ x ')e'i{x') - e i( x ') s "( x ')} - ^ x Wi( x ) e 'i( x ) - e i( x ) s 'i'( x )]^ , (33) 

where the "bulk" pressure has been subtracted, and 

Di = - ^(x'y^x)] 2 - e^iW^x) + S ;(x')e ; (a;)] 2 , (34) 

with the parameter £ being 



t= , (35) 

and J — > is the temporal regulator. This result is obtained either by computing 
the radial-radial component of the stress tensor, or from the total energy. 

In general, this result is divergent. However, consider the special case y/ejl = 
y/e'fi', that is, when the speed of light is the same in both media. Then x = x' and 
the Casimir energy reduces to 

1 f°° °° d 

E = -—j d 2 /e^^(2/ + l)x-ln[l-e 2 (( Si e ; )') 2 ], (36) 



2=1 



where 



If £ = 1 we recover the case of a perfectly conducting spherical shell, treated above. 
In fact E is finite for all £. 

Of particular interest is the dilute limit, wher^3l 

62ira la 
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It is remarkable that the value for a spherical conducting shell, for which £ = 1, is 
only 7% lower, which as Klich remarks}^ is accounted for nearly entirely by the 
next term in the small £ expansion. 

There is another dilute limit which is also quite surprising. For a purely dielectric 
sphere (// = 1) the leading term in an expansion in powers of e — 1 is finite!^^ 

_(._,,>!£!!«!!. ( 39) 

15oo7r a a 
This result coincides with the sum of van der Waals energies of the material making 
up the ball as Ng and I showed earlier in 1998.^3 The term of order (e — l) 3 is 
divergent.^ The establishment of this result was the death knell for the Casimir 
energy explanation of sonoluminescence!^ 



3. Dielectric Cylinders 

The fundamental difficulty in cylindrical geometries is that there is in general no 
decoupling between TE and TM modes. Progress in understanding has therefore 
been much slower in this regime. It was only in 1981 that it was found that the 
electromagnetic Casimir energ y of a perfectly conducting cylinder was attractive, 
the energy per unit length being^lf emjCy i = — 0.01356/a 2 , for a circular cylinder of 
radius a. The corresponding result for a scalar field satisfying Dirichlet boundary 
conditions of the cylinder is repulsivej^ £ n,cyl = 0.000606/a 2 . These ideal limits 
are finite, but, as with the spherical geometry, less ideal configurations have un- 
removable divergences. For example, a cylindrical <5-shell potential has divergences 
(in third order) M^l And it is expected that a dielectric cylinder will have a diver- 
gent Casimir energy, although the coefficient of (s — l) 2 will be finite for a dilute 
dielectric cylinder!^ corresponding to a finite van der Waals energy between the 
molecules that make up the material. In fact, a calculation of the renormalized van 
der Waals energy for a dilute dielectric cylinder gives zero, as is the Casimir en- 
ergy for a cylinder for which the speed of light is the same inside and out, because 
efjb = iP^Ia calculation of the Casimir energy for a dielectric cylinder is in progress 
with my graduate student Ines Cavero-Pelaez. 



4. Conclusions 

We began by reviewing ancient results for the contributions of zero-point energies 
in hadronic (bag) models. We saw that in general there are surface divergences 
in the condensates, and in the energy densities, for the quark and gluon fields. 
We discussed 5- function potentials, which in general give divergent results in 3rd 
order, but are finite in strong coupling. These divergences are largely identified with 
surface energies, which can be interpreted as bulk energies when the boundaries are 
smoothed. Precisely analogous phenomena happen for dielectric balls and cylinders 
(although there is some remarkable symmetry buried in the latter). I hope this 
improved understanding of divergences can lead to improved hadronic models. For 
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the current status of Casimir phenomena, see my review article where complete 
references can be found. 
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